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ABSTRACT 

The fate of ionizing radiation is vital for understanding cosmic ionization, energy budgets in the 
interstellar and intergalactic medium, and star formation rate indicators. The low observed escape 
fractions of ionizing radiation have not been adequately explained, and there is evidence that some 
starbursts have high escape fractions. We examine the spectral energy distributions of a sample of 
local star-forming galaxies, containing thirteen local starburst galaxies and ten of their ordinary star- 
forming counterparts, to determine if there exist significant differences in the fate of ionizing radiation 
in these galaxies. We find that the galaxy-to-galaxy variations in the SEDs is much larger than any 
systematic differences between starbursts and non-starbursts. For example, we find no significant 
differences in the total absorption of ionizing radiation by dust, traced by the 24/mi, 70/im, and 
160/im MIPS bands of the Spitzer Space Telescope, although the dust in starburst galaxies appears 
to be hotter than that of non-starburst galaxies. We also observe no excess ultraviolet flux in the 
GAL EX bands that could indicate a high escape fraction of ionizing photons in starburst galaxies. 
The small Ha fractions of the diffuse, warm ionized medium in starburst galaxies are apparently due to 
temporarily boosted Ho luminosity within the star-forming regions themselves, with an independent, 
constant WIM luminosity. This independence of the WIM and starburst luminosities contrasts with 
WIM behavior in non-starburst galaxies and underscores our poor understanding of radiation transfer 
in both ordinary and starburst galaxies. 

Subject headings: galaxies: evolution - galaxies: ISM - galaxies: starburst - ISM: general - diffuse 
radiation - radiative transfer 



1. INTRODUCTION 

One of the most urgent needs in understanding the 
evolution of cosmic structure is to clarify the origin and 
properties of the background ionizing radiation field. In 
particular, quasar absorption-line systems are now a fun- 
damental constraint on the development of large-scale 
structure in the cosmic web. The statistical proper- 
ties of the Lyman a absorbers appear to be plausibly 
reproduced by numerical simulations of a ACDM uni- 
verse evolving from primordial fluctuations s uch as those 
in the cosmic microwave b ackground (e.g. iDave et"aLl 
Il999t iBi fc Davidsen | [L997L However, uncertainties in 
the background ionizing radiation field are a great ob- 
stacle in exploiting the wealth of absorber data (e.g . 
Lidz et al.|[2006llSokasian et al.ll200l iMadau et al.lll999t 
Rauch Ill998p t seems clear that at redshifts z < 3, the 



UV b ackground is mainly due to QSOs (e.g. lMadau et al.l 
1999). However , at higher redsh i fts, the reduced QSO 
number density (|Fan et al.l 120011: iRichards et"aT1 l2006f) 
implies that another source, presumably s tellar UV ra- 
diation from galaxies (|Bouwens et al.ll2007l ). is the dom- 
inant contributor. The discovery of the large population 
of Lyman-break ga laxies (LBGs) at these redshifts (e.g. 
iSteidel et al.lll996h supports this scenario. 

T here is some evide nce from both lo cal starbursts 
(e.g. lHoopes et al.ll2007ft and in LBGs (e.g. lShaplev et al.l 
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2006) that at least some starbursts have high escape frac- 
tions /esc for ionizing radiation. However, most studies 
find that the / Psr is only a few percent, both locally (e.g. 
i Heckman et al .1120011 : iLeitherer et al. || 1995HGrimes et all 
20071) and at high redshift (e.g iWvithe et al.l 120101: 
ISteidel et all 120011: iShaplev et all 12006ft . iSiana et al l 
(12010ft find that no more than 8% of galaxies at z ~ 1.3 
can have relative escape fractions greater than 0.50, and 
that the average ionizing emissivity appears to decrease, 
approaching z ~ 0. The low observed escape fractions 
are puzzling in vi ew of predictions for h igher values of 
/ csc . For example. iClarke fc Oev I (|2002h predicted that 
above a threshold star-formation rate, the interstellar 
medium (ISM) is shredded by the superwind mechan- 
ical feedback, thereby opening avenues for the escape 
of Lym an continuum radiation. Furthermore, lOev et all 
(2007) found that the fraction of Ho luminosity /wim 
contributed by the diffuse, warm ionized medium (WIM), 
is much lower in starburst galaxies, defined as those 
galaxies with high Ha surface brightnesses. While /wim 
generally has a robust, universal value ~ 0.5 in star- 
forming galaxies, it is systema t ically lower in starbursts, 
by factors of 2 - 5; lOev et a l. (2007) suggested that the 
lower /wim could be evidence of ISM density-bounding 
in the starbursts, suggesting the escape of ionizing radi- 
ation. 

In the present study, we examine the multiwavelength 
spectral energy distributions of star-forming galaxies to 
better understand the fate of the ionizing radiation in 
starbursts vs non-starburst galaxies. Observations in 
three mid- and far-infrared Spitzer wavelengths, the 
MIPS 24/tm, 70/im, and 160/im bands, clarify the role 
of dust heating, while Ha and GALEX NUV and FUV 
bands yield direct estimates of the ionizing stellar popu- 



2 



Hanish et al. 



lation. 

2. DATA 

For our sample, we selected 23 local star-for ming galax- 
ies fro m the sample of 109 galaxies studied by lOev et al.l 
(2007), which correspond to the first data release of the 
Survey for Ionization in Neutral Gas G alaxies (SINGG; 
iMeurer et all 120061: lHanish et all 12006ft . Since a strong 
differentiation was se en in /wtm betwe en the starbursts 
and non-starbursts bv lOev et al.l ((20071 ). we chose 13 rep- 
resentative starburst galaxies and 10 non-starbursts from 
that work for further multi-wavelength study here. We 
note tha t there are severa l definitions of starburst galax- 
ies (e.g. iLee et al.ll2009al) ; we adopt the star formatio n 
intensity (SFI) criterion advocated by iHeckman I ([2005). 
Thus, Ho surface brightness is the corresponding observ- 
able quantity for the SFI, and starbrusts are defined as 
galaxies having Ha surface brightness within the H«- 
derived half-light radius above a threshold log(EH a [erg 
s _1 kpc -2 ]) > 39.4. Similarly, the 10 non-starburst 
galaxies have £hq below this value. 

We obtained new observations and archival data from 
the Spitzer Space Telescope. These observations were 
taken with the MIPS camera, observing in the 24/im, 
70/im, and 160/im channels to estimate the thermal emis- 
sion from hot and cold dust. These MIPS bands are 
diffraction-limited, with effective resolutions of 6, 18, and 
40 arcseconds, respectively. The mosaic images have 
pixel scales of 2.45, 4.0, and 8.0 arcseconds per pixel, 
respectively, although this resolution is the result of a re- 
sampling algorithm applied to images with native pixel 
scales of up to 18 arcseconds per pixel for the 160/im 
source images. As a result, sources in the 70/im and 
160/im images are significantly less resolved than those 
in our other bands. This point spread function (PSF) re- 
sults in a slight underestimation of total infrared fluxes 
due to aperture effects. Aperture correction factors for 
the MIPS bands were estimated bv lDale et all (<2009D to 
range from 1 - 3%, so this effect will not significantly 
alter our observed trends. 

Of the 23 galaxies in our sample, we obtained new 
Spitzer observational data for eight, as shown in Ta- 
ble [TJ Twelve additional galaxies had been fully observed 
by previous studies, four of which had archival data sets 
from multiple previous sources. As a result, twenty of 
our sources (11 starbursts, 9 non-starburst galaxies) pos- 
sessed the full array of MIPS observations, either from 
archival data or from our own requested observations. 
The three galaxies with incomplete MIPS data were NGC 
178 (lacked 24/im and 160/im data), NGC 1808 (lacked 
160/im data), and NGC 3365, which had no MIPS data 
at all due to the termination of the cold Spitzer mission. 

We used available Ha and i?-band data from SINGG 
Release 1. This survey uses an R band continuum filter 
centered at 6508 A and a narrow-band Ha filter appropri- 
ate to the redshift of each galaxy to create a continuum- 
subtracted emission-line image. The result is an accurate 
measure of Ha line flux, tracing the formation of high- 
mass stars within each galaxy. SINGG images have a 
pixel resolution of 0.432 arcseconds per pixel, although as 
a ground-based survey, the images are limited by the ef- 
fects of atmospheric seeing, which generally ranged from 
1.0 to 1.5 arcseconds. The SINGG field is approximately 
14 arcminutes on each side; only one of our sources, NGC 



5236 (Messier 83), was found to contain features extend- 
ing beyond the bounds of its SINGG images. 

We acquired archival Galaxy Evolution Explorer 
(GALEX) data for all 23 galaxies. These data consisted 
of two UV bands, 1516 A FUV (covering a range from 
1350 A to 1750 A) and 2271 A NUV (ranging from 1750 
A to 2800 A); while actual ionizing radiation remains 
shortward of both bands, these serve as approximate 
indicators of ionizing flux. The detector used for both 
bands has a pixel scale of 1.5 arcseconds per pixel and a 
PSF FWHM of 4.5 - 6 arcseconds, resulting in less de- 
tail than the SINGG images, but significantly more than 
the Spitzer bands. GALEX images cover a far larger 
field than any of our other bands, a circular field with a 
diameter of approximately 1.2 degrees; as a result, the 
GALEX NUV image served as the positional reference 
for the alignment algorithms used in our masking scripts. 

GALEX and SINGG data required no additional data 
reduction. However, the Spitzer observations required 
reprocessing with different parameters to ensure a bet- 
ter fit of the sky background for extended sources. The 
method used to generate the intermediate Basic Cali- 
brated Data (BCD) images, known as the Germanium 
Reprocessing Tools (GeRT), defaults to rejecting a 5x5- 
pixel box around the center point of each object before 
calculating a sk y profile. We instead u sed Source Ex- 
tractor (v. 2. 5.0, IBertin fe Arnouts"l lT996) to identify the 
visibly extended features of each galaxy; this allowed us 
to then mask a larger area around each of these regions 
with a circular mask whose size is based on the semimi- 
nor axial length of each galaxy, resulting in a significantly 
improved model fit to the sky background. Once this was 
done, the MOPEX software package was used to combine 
the modified BCD images into a single mosaic for each 
galaxy; only minor alterations to the default MOPEX 
process were needed, such as compensating for differing 
sky background levels in those few cases where we com- 
bined MIPS data from different archival sources. 

All data processing beyond this stage was performed 
with IDL, the Interactive Data Language; all scripts 
were either custom-mad e or dr awn from the ID L As - 
tronomy User's Library (http://idlastro.gsfc.nasa.gov/), 
more commonly known as "astro-lib" . For each galaxy, 
the various bands were aligned using IDL scripts con- 
tained in astro-lib. Foreground objects were removed 
with a common mask; an initial mask was derived by 
applying Source Extractor to the Spitzer 24/im and 
GALEX NUV images, and then iteratively modified by 
user input based primarily on a 3-color image similar to 
that of Figure QJa), combining images from all three tele- 
scopes, while a series of alternate 3-color images (such 
as those telescope-specific images shown in Figure QIb) 
through (d)) were used to evaluate any questionable ob- 
jects. The resulting mask was then scaled and rotated 
appropriately to match each input image's coordinates. 
The Spitzer 70/im and 160/im images were not signif- 
icantly contaminated by foreground objects, as stellar 
continua are extremely faint at far-infrared wavelengths, 
so we did not use masks in these bands. The end result 
was eight bands of data, each containing a pixel mask for 
foreground objects derived from the same source, allow- 
ing for a consistent flux measurement process. 

The flux within each band was measured within a sin- 



Fate of Ionizing Radiation in Starbursts 



3 




Figure 1. 3-color images, using NGC 1421 (top) and NGC 1291 (bottom) as examples of our starburst and non-starburst galaxies, 
respectively. Columns shown are: (a) The primary composite image used in our masking: Spitzer 24/mi flux is red, SINGG R flux is green, 
and GALEX NUV is blue; (b) the GALEX bands: NUV flux is red, FUV flux is blue, and the average of the two bands is green; (c) the 
SINGG bands: R flux is blue, narrow-band flux is green, and the continuum-subtracted Ha flux is red; (d) the Spitzer bands: 24/^m flux 
is blue, 70fj,m flux is green, and 160^tm flux is red. Images of this type for all 23 galaxies are available in the online edition of this paper. 



gle elliptical aperture for each galaxy, whose center posi- 
tion, position angle, and axial ratio were taken directly 
from the SINGG project. The fluxes were measured us- 
ing the original input images, instead of the rotated and 
scaled versions used to generate the pixel masks. How- 
ever, the SINGG apertures were derived based on R band 
flux and extended Ha emission, and so they tended to 
be large relative to the measurable size of the Spitzer 
band emissions. To maintain consistency, we limited our 
measurement of flux ratios to that part of each galaxy 
within the half-light radius measured in Ha, r^o, and 
this extended to derived quantities such as Ha equiva- 
lent width. 

Corrections of internal extinction were drawn from 
relationships involving total infrared flux to Ha flux 
or UV flux, as appropriate; the GALEX extinction 
prescriptions. Eq uations Q] and [2j were drawn from 
iBuat ct al. (2005). T he ex tinction relation for Ha was 
from lKennicutt et al.l ([2009) , shown in Equation^ while 
the i?-band extinction was assum ed to be directl y pro- 
portional to the Ha extinction by iCalzetti I ([2001). The 
extinction relations, in magnitudes, were 

x = log(Fm/ (f/^NUv)) 

A NVV = -0.0495x 3 + 0.4718a; 2 + 0.8998x + 0.2269 (1) 

J/ = log(Jk./(i// v> Fuv)) 

Afuv = -0.0333j/ 3 + 0.3522y 2 + 1.1960y + 0.4967 (2) 

A Ha = 2.5 log(1.0 + 0.0024F ir /(F Hq )) (3) 

Aa = A H a/2.3 (4) 

As our infrared fluxes primarily represent thermal 
emission from various temperatures of gas and dust, no 
extinction corrections were applied to our three MIPS 
bands. All of the above equations link dust extinction 
in each band to the total infrared flux. To estimate the 
total Spitzer-derrved dust emission, we calcul ated an in- 
tegrated infrared flux from the model fits of iDale et al.l 



(2007), using their derived relationship 

Fir = 1.559 !//,,24 M m +0.7686 vf lr ,7o ltm + 1.347 vf v , 160/im (5) 

This expression, representing the total infrared flux be- 
tween 3/zm and llOO^m, corresponds primarily to the 
thermal emission of dust within each galaxy. All fluxes 
used to determine extinction corrections were measured 
using the full SINGG-defined apertures instead of the 
reduced Ha-derived apertures used for all other fluxes 
in this paper, to maintain consistency with the studies 
which derived these relationships. 

The extinction-corrected fluxes are given in Table [2j 
column 1 lists the galaxy ID, and columns 2-8 give 
our measured fluxes in the seven bands we studied. The 
Ha flux densities are calculated from the emission line 
fluxes by assuming a constant line width of 1.026 A, cor- 
responding to the thermal broadening of a source at a 
temperature of 12,000 K. While fluxes were directly avail- 
able for several of these sources from the projects which 
provided many of the data images, we measured all fluxes 
ourselves to ensure consistent masks, apertures, and in- 
tegration techniques were used. The listed uncertainties 
are those due to sky measurement. While many other 
sources of error exist, such as those of the [Nil] correction 
and continuum scaling ratio used in the derivation of Ha 
line flux, the total uncertainties are generally domi nated 
by the extinction model used ([Hanish et al.l 120061 ). As 
Equations |TJ — 0] have no quoted uncertainties, we cannot 
easily quantify this effect. 

3. COMPARISON OF STARBURST AND NON-STARBURST 
GALAXIES 

We calculate the star formation rates for each galaxy 
by three methods. Our UV-d erived and opticall y-derived 
SFRs use the relationships of Kcnnicutt (1998), and our 
infrared SFR uses that of ICalzetti et al.l (|2007l) . These 
relationships are, for a Salpeter IMF, 

SFRuv[A4 yr- 1 ] = 1.4 x 10" 28 L„, FU v[erg s^Hz" 1 ] (6) 
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Table 1 

Aperture definitions, by galaxy 



Galaxy 



RA (J2000) 
[h:m:s] 



Dec (J2000) 
[d:m:s] 



Distance 
[Mpc] 



log(-M H i) 
[M Q ] 



Vhel 

[km s _ 



50 (Ha) 



PA 
[deg] 



a/b 



Spitzer 
reference^ 



GALEX 



survey 



Starburst Galaxies 



ESO409-IG015 


0:05:31.7 


-28:05:49.2 


10.24 


8.27 


737±6 


11.5 


141.1 


1.72 


(1) 


(7) 


NGC178 


0:39:08.2 


-14:10:26.4 


19.8 


9.40 


1447±3 


20.8 


8.6 


2.12 


(1) 


(7) 


NGC625 


1:35:03.1 


-41:26:13.2 


4.45 


8.09 


396±1 


43.9 


94.0 


2.96 


(2) 


(7) 


NGC922 


2:25:03.8 


-24:47:27.6 


41.2 


10.07 


3082±5 


30.7 


176.0 


1.09 


(1) 


(7) 


NGC1421 


3:42:29.3 


-13:29:20.4 


27.7 


9.85 


2087±5 


55.2 


177.4 


3.26 


(1) 


(8) 


NGC1487 


3:55:45.6 


-42:22:01.2 


10.1 


9.25 


848±1 


31.3 


60.5 


2.52 


(2) 


(7) 


NGC1510 


4:03:32.6 


-43:24:00.0 


11.0 


9.54 


913±10 


3.3 


129.4 


1.13 


(2,3) 


(9) 


NGC1705 


4:54:14.2 


-53:21:39.6 


6.46 


7.96 


633±6 


20.7 


44.1 


1.38 


(3,4) 


(9) 


NGC1800 


5:06:25.2 


-31:57:18.0 


9.7 


8.54 


807±1 


16.1 


114.1 


1.52 


(2) 


(9) 


NGC1808 


5:07:43.2 


-37:30:32.4 


12.3 


9.53 


995±4 


39.0 


134.3 


2.10 


(5) 


(7) 


NGC5236 


13:37:00.0 


-29:51:56.2 


4.03 


9.89 


513±2 


122.0 


91.3 


1.02 


((>) 


(9) 


NGC5253 


13:39:55.4 


-31:38:27.6 


2.54 


8.22 


407±3 


12.6 


44.8 


2.01 


(2) 


(9) 


NGC7126 


21:49:18.2 


-60:36:36.0 


42.7 


10.46 


2981±2 


17.8 


67.4 


1.73 


(1) 


(7) 


N on-Starbur 


st Galaxies 




















NGC1291 


3:17:18.0 


-41:06:32.4 


10.1 


9.41 


839±2 


133.8 


76.7 


1.17 


(3) 


(9) 


NGC1311 


3:20:06.7 


-52:11:06.0 


5.95 


8.25 


568±5 


31.8 


38.6 


2.60 


(2) 


(10) 


IC2000 


3:49:07.4 


-48:51:32.4 


11.9 


9.10 


981 ±4 


54.5 


81.9 


2.97 


(1) 


(10) 


NGC1512 


4:03:54.2 


-43:20:52.8 


10.8 


9.64 


898±3 


69.6 


52.8 


1.04 


(2,3) 


(9) 


NGC1744 


4:59:58.1 


-26:01:22.8 


8.94 


9.56 


741±2 


89.0 


164.6 


1.68 


(2) 


(11) 


ESO486-G021 


5:03:19.7 


-25:25:26.4 


10.3 


8.60 


835±3 


23.1 


97.0 


1.64 


(2) 


(8) 


UGCA175 


9:43:36.5 


-5:54:43.2 


30.8 


9.36 


2026±5 


35.7 


20.4 


1.29 


(1) 


(7) 


NGC3365 


10:46:12.7 


+1:48:50.4 


13.3 


9.18 


986±1 


45.8 


157.2 


3.25 


N/A 


(7) 


UGCA307 


12:53:56.6 


-12:06:14.4 


8.0 


8.67 


821±1 


21.7 


142.6 


2.14 


(1) 


(7) 


UGCA320 


13:03:17.3 


-17:25:26.4 


7.04 


9.12 


742±2 


42.0 


117.4 


3.36 


(2) 


(7) 



Note. — f Refer ences: (1) Our request ed data, Spitzer GO-50332; (2) Spitzer GO-40204 (|Kennicutt et al. 2007); 
(3) Spitzer GO-159 l|Kennicutt et al.1 [20031 ): (4) Spitzer GO-20528 (PI: Martin, C); (5) Spitzer GO-86 (PI: Werner, 
M.); (6) Spitzer GO-59 (PI: Rieke, G.) 

J Surveys: (7) Guest Investigators Survey; (8) All-sky Imaging Survey; (9) Nearby Galaxy Survey; (10) Medium Imaging 
Survey; (11) LGAL Survey 



Table 2 

Flux densities (log(/„ [Jy])^) within r5o(Ha), by galaxy 



Galaxy 




FUV 






NUV 






R 






Ha 






24/im 






70/^m 






160/im 


Starburst Galaxies 








































ESO409-IG015 


-3 


008±0. 


.004 


-3 


.076±0 


.004 


-2 


,783±0. 


.003 


-0 


489±0. 


.001 


-1 


.923±0. 


007 


-1 


,191±0. 


012 


-1 


.603±0 


.017 


NGC178 


-2 


456±0. 


.001 


-2 


.413±0 


.001 


-1 


,780±0. 


.001 


-0 


,214±0. 


.001 




N/A 




-0 


194±0. 


.006 




N/A 




NGC625 


-1. 


,849±0. 


001 


-1 


.881±0 


.001 


-1 


230±0. 


.002 





,581±0. 


.001 


-0 


.526±0. 


004 





290±0. 


.009 





.202±0 


Oil 


NGC922 


-1 


,822±0. 


001 


-1 


819±0. 


.003 


-1 


364±0. 


.003 





398±0. 


001 


-0 


.543±0. 


003 





530±0. 


.005 





.417±0 


.005 


NGC1421 


-1 


831±0. 


002 


-1 


,780±0. 


.003 


-1 


032±0. 


001 





605±0. 


003 


-0 


,246±0. 


002 





,810±0. 


.005 





.854±0 


.004 


NGC1487 


-1. 


953±0. 


001 


-1 


959±0. 


.001 


-1. 


440±0. 


002 





.329±0. 


001 


-0 


,856±0. 


004 


0. 


164±0. 


006 





.044±0 


.010 


NGC1510 


-2 


,758±0. 


001 


-2 


862±0. 


.001 


-2 


169±0. 


001 


-0 


231±0. 


001 


-1 


409±0. 


001 


-0 


980±0. 


.006 


-1 


.323±0 


.015 


NGC1705 


-1. 


,849±0. 


001 


-1 


,885±0. 


.001 


-1 


.496±0. 


001 





,275±0. 


001 


-1 


,467±0. 


004 


-0 


.316±0. 


003 


-0 


599±0 


.005 


NGC1800 


-2 


333±0. 


001 


-2 


,357±0. 


.001 


-1. 


696±0. 


002 


-0 


190±0. 


002 


-1 


638±0. 


010 


-0. 


558±0. 


.006 


-0 


.621±0 


Oil 


NGC1808 


-2 


,078±0. 


001 


-1 


,832±0. 


.001 


-0. 


615±0. 


001 


1 


191±0. 


001 





703±0. 


001 


1. 


474±0. 


006 




N/A 




NGC5236 


-0 


,389±0. 


001 


-0 


.385±0 


.001 


0. 


.425±0. 


006 


2 


020±0. 


.010 


1 


.431±0. 


001 


2 


,144±0. 


.006 


2 


.436±0 


.015 


NGC5253 


-1 


,378±0. 


001 


-1 


,594±0. 


.001 


-1 


311±0. 


.003 


1 


,241±0. 


001 





.561±0. 


001 





.512±0. 


.004 





.158±0 


.005 


NGC7126 


-2. 


.561±0. 


001 


-2 


503±0. 


.002 


-1. 


433±0. 


001 





,084±0. 


001 


-0 


,907±0. 


002 





079±0. 


.004 





.024±0 


.004 



NGC1291 

NGC1311 

IC2000 

NGC1512 

NGC1744 

ESO486-G021 

UGCA175 

NGC3365 

UGCA307 

UGCA320 



-2.091±0.004 
-2.626±0.001 
-2.422±0.004 
-1.926±0.002 
-2.030±0.001 
-2.855±0.007 
-2.617±0.006 
-3.150±0.006 
-3.071±0.002 
-2.671±0.001 



-1.844±0.003 
-2.542±0.002 
-2.389±0.006 
-1.878±0.003 
-1.990±0.002 
-2.796±0.014 
-2.579±0.011 
-2.958±0.010 
-3.078±0.004 
-2.602±0.002 



0.171±0.002 
-1.841±0.002 
-1.424±0.002 
-0.525±0.002 
-1.069±0.005 
-2.070±0.007 
-1.697±0.002 
-1.753±0.002 
-2.604±0.007 
-2.040±0.004 



0.607±0.031 
-0.296±0.003 
-0.146±0.012 
0.417±0.020 
0.265±0.084 
-0.584±0.011 
-0.422±0.017 
-0.470±0.007 
-0.579±0.003 
-0.111±0.004 



-0.575±0.083 
-1.907±0.044 
-1.218±0.040 
-0.528±0.033 
-1.098±0.090 
-2.109±0.072 
-1.480±0.032 

N/A 
-2.261±0.047 
-1.899±0.078 



0.541±0.023 
-0.650±0.013 
-0.226±0.011 
0.597±0.014 
0.102±0.023 
-0.896±0.018 
-0.431±0.023 

N/A 
-1.507±0.066 
-0.710±0.026 



0.904±0.018 
-0.735±0.019 
0.129±0.005 
0.954±0.017 
0.530±0.018 
-0.986±0.037 
-0.193±0.008 

N/A 
-1.435±0.033 
-0.892±0.057 



Note. — f All flux densities have been corrected for internal extinction, and quoted uncertainties only include 
sky measurement errors. 



Fate of Ionizing Radiation in Starbursts 



5 



SFR H c<[-Ma yr' 1 ] = 7.9 x Kr 42 L Ha [erg s" 1 ] (7) 
SFR m [X yr- 1 ] = 2.02 x 1(T 38 K„, 2 4 Mm [erg s' 1 ]] a88 f8) 

Quantities derived from our directly measured values 
are given in Tabled Column 1 lists the galaxy ID, and 
column 2 shows the fraction of Ho luminosity /wim ob- 
served in t he diffuse, warm ionized medium, taken di- 
rectly from IQev et al.l (|2007l ). Column 3 gives the to- 
tal infrared luminosity for each galaxy based on Equa- 
tion [3 column 4 gives the equivalent width of the Ha 
emission line within the Ha effective radius, column 5 
shows the ratio of the UV-derived SFR from Equation [6] 
to the Ha-derived SFR from Equation [Jj and column 
6 gives the concentration index of our galaxies in Ha, 
defined as the ratio of the semimajor axis of the aper- 
ture containing 50% of the Ha flux to that of the aper- 
ture containing 90%. Appropriate extinction corrections 
have been applied to the values of equivalent width (as 
this is just Pka/ Jv.r) and log(SFR(Ha)/SFR(UV)) given 
in columns 4 and 5. The mean and standard deviation 
of each subsample is given for each variable, illustrating 
the differences between starburst and ordinary galaxies 
in quantities like /wim and Ha equivalent width. 

All of the extinction corrections in Equations Q] - 0] in- 
volve the ratio of Fir to the flux at the appropriate wave- 
lengths; while the IR luminosities given in Table|3]tend to 
be somewhat higher for starburst galaxies than for non- 
starbursts, the fluxes at other wavelengths increase by 
comparable amounts, resulting in similar distributions 
of extinction corrections. For instance, the Ha extinc- 
tion corrections of Equation [3] have a mean of 0.49 ± 0.32 
magnitudes for starburst galaxies and 0.38 ± 0.23 mag- 
nitudes for non-starbursts, although this is the result 
of very inhomogenous populations; for instance, four 
out of thirteen starburst galaxies have Ha corrections 
above 0.95 magnitudes, but the remaining nine all pos- 
sess corrections of 0.45 magnitudes or less. The Ha cor- 
rections for ordinary star-forming galaxies range from 
0.09 - 0.71 magnitudes. The R band extinctions are, 
as given in Equation [4j directly proportional to the Ha 
extinctions and so follow a similar pattern. The ultra- 
violet extinctions tend to be slightly higher than those 
of Ha, especially for starburst galaxies; in magnitudes, 
the NUV extinctions were 0.73 ± 0.53 for starbursts and 
0.41 ± 0.40 for ordinary galaxies, while the FUV extinc- 
tions were 1.01 ± 0.67 for starbursts and 0.68 ± 0.58 for 
non-starbursts. 

3.1. Spectral Energy Distributions 

Figure [2] shows the spectral energy distribution (SEDs) 
for the galaxies in our sample, normalized by the Ha line 
flux, an indicator of star formation rate. Starburst galax- 
ies are shown with solid lines, and non-starbursts with 
dashed. The mean composite SEDs for the starbursts 
and non-starbursts are shown with the heavier black solid 
and dashed lines, respectively, with error bars displaying 
the standard deviations of mean at each wavelength. The 
variation between the individual galaxies is much larger 
than any systematic variation between the starburst and 
ordinary galaxy subsets. That said, the starburst galax- 
ies produced an average of 0.07 dex less ultraviolet flux 
per unit of Ha than their ordinary counterparts. Like- 
wise, for the 20 galaxies with complete MIPS data, the 



starburst galaxies produced an average of 0.03 dex less 
total infrared flux per unit Ha. While at 24/im we mea- 
sured a difference of almost 0.3 dex between starbursts 
and ordinary galaxies, this gap was effectively negligible 
at 70/^m, and at 160/im there was an anticorrelation of 
over 0.3 dex, resulting in a net discrepancy of 0.03 dex 
in total infrared flux from Equation [5j with starbursts 
again producing less flux per unit of Ha than their or- 
dinary counterparts. This difference is far less than the 
measurement uncertainties in each point, and is negligi- 
ble compared to the statistical variation between the data 
points. In both IR and UV, the individual galaxies' flux 
ratios often varied by up to 0.7 dex from the mean. We 
note, however, that some of this scatter is introduced by 
galaxies with low Ha fluxes, which tend to have flux ra- 
tios further from the mean and significantly higher frac- 
tional uncertainties from sky subtraction. Within our 
sample, these galaxies tended to have infrared-to-Ha ra- 
tios well below those of their larger counterparts, imply- 
ing that our mean flux ratios might be slight underesti- 
mates. 

Figure [3] shows a modest correlation between the 
24/im/70/im flux ratio and the Ha/UV star formation 
rate ratio. Instead of showing the Ha/UV flux ratio 
itself, the abscissa in Figure [3] shows the ratio of star- 
formation rates derived using Equations [5] and which 
are proportional to the SINGG Ha and GALEX FUV 
fluxes respectively. As noted above, Figure [2] does show 
that the starbursts tend to have larger 24/!m/70/im flux 
ratios; as this ratio should be expected to increase with 
temperature for purely thermal emission, this supports 
the observ ation that starburs ts possess hotter dust tem- 
peratures (jDraine fc Li Il2007t) , although this discrepancy 
would be far less without the presence of NGC 5253, the 
starburst galaxy with by far the highest 24/im/70/im ra- 
tio in our sample. A correlation with the Ha/UV ratio 
links the 24/im/70/im ratio directly to the starburst ac- 
tivity itself. 

The Ha/UV ratio can, in theory, be used to check 
whether ionizing photons are preferentially leaking into 
the IGM. In practice, however, such an analysis is made 
difficult by the fact that many other processes can lead 
to an Ha/UV ratio that is lower than expected. A 
number of studies based upon recent GALEX data (e.g . 
iMeurer et alj 120091 : iLee et al.ll2009ri iBoselli et ai1l2009i) 
as w e ll as on earlier UV datasets (|Bell fc K ennicutt 
2001; Sullivan et al.l [20041 ) have examined this issue in 
detail, and considered effects ranging from variations 
in the short-timescale star formation histories and dust 
attenuation, to the escape of ionizing photons and an 
IMF which is deficient in the most massive stars. It 
is likely that most of the scatter of the Ha/UV ra- 
tio in Figure [3] is due to non-constant star formation 
rates. Nevertheless, it is notable that the average SFR 
ratios for the starburst and non-starburst samples are 
similar (< log(SFR(Ha)/SFR(UV))> = 0.06 ±0.30 and 
0.02 ±0.18, respectively, as shown in Table[3]), with sim- 
ilar ranges of values. This argues that processes which 
lead to variations in the Ha/UV ratio, inc luding Lyman 
continuum photon escape (|Lee et al.ll2009bf ). do not seem 
to be systematically dominant in the starburst popula- 
tion relative to ordinary star-forming galaxies. 
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Table 3 

Derived flux values, by galaxy 



Galaxy 


/wim' 1 ' 


log(L TIR [crg s x ]) 


log(EW Ha [A]) 


locf SFR( h 

iL *&V SFR CTTV1 / 


7^( Hq ) 


Starburst Galaxies 










ESO409-IG015 


0.15 


40.80 


2.29 


0.13 


0.67 


NGC178 


0.37 


N/A 


1.57 


-0.15 


0.41 


NGC625 


0.21 


41. oy 


i si 
l.ol 


U.U4 


Ki.OI 


NGC922 


0.29 


43.68 


1.76 


-0.17 


0.77 


NGC1421 


0.47 


43.67 


1.64 


0.04 


0.70 


NGC1487 


0.40 


42.10 


1.77 


-0.11 


0.35 


NGC1510 


0.22 


41.25 


1.94 


0.13 


0.47 


NGC1705 


0.52 


41.16 


1.77 


-0.27 


0.35 


NGC1800 


0.50 


41.35 


1.51 


-0.25 


0.33 


NGC1808 


0.19 


N/A 


1.81 


0.87 


0.34 


NGC5236 


0.31 


43.52 


1.59 


0.02 


0.58 


NGC5253 


0.11 


42.22 


2.55 


0.23 


0.30 


NGC7126 


0.54 


43.30 


1.52 


0.25 


0.24 


Mean 2 


0.33±0.15 


42.2±1.1 


1.81±0.31 


0.06±0.30 


0.47±0.17 


Non-Starburst Galaxies 










NGC1291 


0.91 


42.66 


0.44 


0.30 


0.42 


NGC1311 


0.65 


40.79 


1.55 


-0.06 


0.35 


IC2000 


0.80 


42.05 


1.28 


-0.12 


0.50 


NGC1512 


0.72 


42.77 


0.94 


-0.05 


0.42 


NGC1744 


0.61 


42.14 


1.33 


-0.10 


0.40 


ESO486-G021 


0.70 


41.02 


1.49 


-0.12 


0.56 


UGCA175 


0.66 


42.60 


1.28 


-0.20 


0.47 


NGC3365 


0.71 


N/A 


1.28 


0.29 


0.38 


UGCA307 


0.59 


40.37 


2.03 


0.10 


0.20 


UGCA320 


0.57 


40.86 


1.93 


0.17 


0.24 


Mean 2 


0.69±0.10 


41.7±0.9 


1.35±0.46 


0.02±0.18 


0.39±0.11 



Note. — (1) From [Oev et al.l ((2007); (2) Mean and standard deviation of each sub- 
sample. 



160um 70um 24um Ha NUV FUV 




12.5 13.0 13.5 14.0 14.5 15.0 

log(v [Hz]) 

Figure 2. The SEDs for the galaxies in this sample. Starburst galaxies use solid lines, while ordinary galaxies are represented by dashed 
lines, with line colors reflecting each galaxy's total IR flux ranging from high (red) to low (blue). The heavier black lines are the averages 
for these two types. 



3.2. Radial Profiles 



We also attempted to determine the mean radial dis- 
tribution of the flux within each band, in order to com- 
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log(SFR(Ha) / SFR(UV)) 
Figure 3. Relationship between the ratio of 24fim flux density 
to that at 70/xm and the ratio of Ha star formation rate to that 
measured using UV flux. Circles denote starbursts, diamonds are 
ordinary galaxies. 

pare the concentrations of star-forming regions to those 
of the hot and cold dust. To do this, we deconvolved 
each IR and UV image by an appropriate PSF through a 
maximum entropy algorithm; this had the greatest effect 
on the Spitzer MIPS images, due to their poor spatial 
resolutions. We then measured the flux distribution of 
the resulting deconvolved image using the same annular 
method as in the earlier measurement of the total flux. 
The results were plotted in Figure [4] as the percentage 
of the total flux F max as a function of the percentage of 
galaxy radius r max , for each band. The composite star- 
burst and non-starburst radial profiles are shown by solid 
and dashed lines, respectively, with the values of r rnax 
determined from the SINGG Ha data for each galaxy. 

The resulting profiles show several interesting trends. 
The two ultraviolet bands are nearly identical in their ra- 
dial distributions for both starbursts and ordinary galax- 
ies, and both match well to the optical i?-band distri- 
butions. The link between Ha and 24/im emission in 
starburst galaxies, suggested by Figure [31 is further con- 
firmed by the coincidence of their radial distributions. 
Figure H] shows that in starbursts, Ha and 24/im emis- 
sion are both strongly concentrated, with similar radial 
profiles, suggesting that the hot, young stellar popula- 
tion is directly responsible for the warmer dust traced 
by the 24^tm emission. Similarly the i?-band continuum 
correlates well to the UV bands for both starbursts and 
non-starbursts, with both acting as tracers of stellar pop- 
ulation. 

On the other hand, 24/im flux for non-starburst galax- 
ies is the least concentrated of any band. Both 70/im and 
160/im infrared fluxes were, on average, less concentrated 
than the other bands, but the data in these bands do 
not lend themselves well to measuring accurate profiles 
at small radii due to their low spatial resolutions, small 
fields of view, an imprecise deconvolution method, and 
uneven sky backgrounds. While these factors are not sig- 
nificant enough to change the general trends seen in Fig- 
ure IH they do explain why the MIPS contours are more 
uneven than those of the SINGG and GALEX bands. 

Given the observed correlation between the Ha and 
24fim radial profiles, we then evaluated whether the to- 
tal IR flux implies significant absorption of ionizing radi- 
ation by dust in starburst galaxies. While photons may 
be absorbed by dust at a wide range of wavelengths, 
the high-energy photons capable of ionizing hydrogen are 



more likely to be absorbed by even a moderate quantity 
of dus t or gas and are unlikely to escape a galaxy en- 
tirely ([Lee et all l2009bj) . If this hypothesis is correct, 
we would expect the total infrared flux to be larger in 
starburst galaxies relative to their Ha line flux, and the 
starbursts with more intense star formation would tend 
to have lower rso/Vgo ratios, since their star formation 
activity is more concentrated. Figure [5] shows the re- 
lationship between the IR/Ha flux ratio and the radial 
concentration of Ha. We find no significant difference 
between starbursts and non-starbursts in this sample, 
either in the IR-to-Ha ratio or in concentration index; 
the mean log(FiR,/FH Q ) for starbursts was 1.82 ± 0.37 
while our non-starburst galaxies had 1.72 ±0.43, and the 
distributions of r§o /Vgo values are given in Table [3] 

3.3. The Warm, Ionized Medium 

As mentioned in § [TJ IQev et all (|2007|) suggested that 
the small WIM fractions /wim seen in Ha observations of 
starburst galaxies could be caused by the escape of ioniz- 
ing radiation. If the starburst activity yields enough ion- 
izing radiation that all of the diffuse ISM is photoionized, 
then even if the starburst luminosity increases, the total 
WIM luminosity remains constant or decreases, causing 
/wim to decrease. High / csc is furthermore plausible in 
starburst galaxies, based on th eoretical predictions. For 
example, iClarke fc Oev I (|2002l ) suggested the existence 
of a threshold star formation rate above which a galaxy's 
ISM is shredded by mechanical feedback, allowing UV 
radiation to escape freely: 

CT7-D r a a -li n1 r*flSM,lO% /„S 

SFR cni [M yr J = 0.15 ; (9) 

Id 

where Mjsm 10 is the ISM mass in units of 10 10 -M , 
vio is the thermal velocity dispersion of the ISM in units 
of 10 km s _1 , and fd is a geometric correction factor 
for disk galaxies. This relation results from a simple 
criterion that balances the supernova mechanical energy 
resulting from star formation against the total ISM ther- 
mal energy. If the former dominates, the ISM is shred- 
ded, a galactic outflow is generated, and ionizing photons 
escape. As detailed in IQev et al.l ()2007f ). for this discus- 
sion the ISM mass can be approximated by the HiPASS- 
derived Hi mass Aim, the thermal velocity dispersion of 
the ISM is assumed to be a constant 10 km s _1 . fd is 
assumed to be 0.1 for our entire sample to account for 
a flattened, disk geometry; as a result of this geometry 
the star form ation intens i ty, sh own to relate to galac- 
tic outflow by iHeckmanl ((2002), can be assumed to be 
directly proportional to the total star formation rate of 
each galaxy. Taken together, these assumptions make 
SFR cr i t directly proportional to the Hi mass. 

The SFR derived from each band, according to Equa- 
tions [fO and [51 are shown in Figure O normalized by 
Hi mass. The starburst galaxies are the only ones with 
star formation rates consistently surpassing the SFR cr it 
line given by Equation [9l with nearly all starbursts hav- 
ing at least one tracer of star formation fall within 0.2 
dex of this critical threshold. In contrast, almost all ordi- 
nary star-forming galaxies fall at least 0.3 dex below the 
SFR cr i t line, with only UGCA 175 having even a single 
tracer of star formation with a value above the critical 
threshold. This confirms our earlier findings based only 
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Figure 4. Radial flux profiles in each wavelength band, using deconvolved images. Solid lines are the means of the profiles for starburst 
galaxies, dashed lines represent ordinary star-forming galaxies. The different bands are color-coded as shown. 
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Figure 5. Ratio of integrated total IR flux to Ha line flux from 
Table [21 as a function of concentration index rso/rgo. Circles de- 
note starbursts, diamonds are ordinary galaxies. 

on Ha observations, that the starburst galaxies are ex- 
pected to have high / csc , although we caution that the 
crude relation in Equation [9] needs refinement. Another 
possibility is that SFR cr it is underestimated if the mass 
of the ISM is predominantly molecular. 

However, in § 13. 1[ we saw that the Ha-derived star- 
formation rates for starburst galaxies are not systemat- 
ically suppressed relative to those determined from UV 
or IR data, as would be expected for high / esc . The 
SFR calibrations are based on total luminosities, rather 
than star-formation intensity, which is the differentiating 
factor here between starbursts and massive star-forming 
galaxies. Thus, the similarity between starburst and non- 
starburst SEDs suggests no significant / esc in the star- 
bursts. 

We therefore seek an alternative explanation for the 
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Figure 6. Plotting SFR/AIhi versus the critical SFR needed to 
shred the ISM. Circles represent IR, diamonds are UV, squares 
are Ha and the dotted line is the Hi-based SFR cr it derived from 
equation [2] Filled symbols represent starbursts, open symbols are 
ordinary galaxies. 
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low /w tm observed in the starburst galaxies of lOev et al.l 
(2007). One possibility is a disproportionate absorp- 
tion of Lyman continuu m photons by dust in the star- 
bursts (|Lee et al.ll2009b[ ). However, Figure [7] shows no 
correlation between the deviation from the mean in the 
IR/Ha flux ratio and the fraction of diffuse, warm ionized 
medium (/wtm), which we know to be lower in starburst 
galaxies. This is also expected from our inference above, 
that any inferred difference between the starbursts and 
non-starbursts is insignificant for this sample. 

A more plausible explanation for the high WIM frac- 
tions in starbursts is suggested by Figure [51 which shows 
/wim as a function of Ha equivalent width: a strong 
anticorrelation is apparent. The Ha equivalent width is 
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Figure 7. Plotting the deviation from the mean of the (IR/Ha) 
ratio for each galaxy in our sample against the fraction of warm, 
ionized medium derived from our earlier work. Circles represent 
starbursts, diamonds are ordinary galaxies. 
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Figure 8. Ha equivalent width as a function of the warm ionized 
medium fraction (/wim)- Circles denote starbursts, diamonds arc 
ordinary galaxies. 

essentially a measure of the starburst age, since the line 
emission is strongest in the youngest regions, while the 
stellar continuum changes on much longer timescales and 
increases with age. Thus, the anticorrelation in Figure [8] 
strongly suggests that starburst galaxies simply have a 
temporarily augmented, intense star formation, which in- 
creases the total Ha luminosity concentrated within the 
Hll regions. In contrast, the luminosity of the diffuse 
WIM remains largely independent of the starburst ac- 
tivity, and so its fractional contribution to a starburst 
galaxy's total Ha luminosity decreases. 

Table [3] shows that, as expected, the mean starburst 
has an Ha equivalent width 0.46 dex higher than that 
of the mean ordinary galaxy, corresponding to a factor 
of three increase in mean star formation intensity. In 
fact, there was very little overlap between the two sam- 
ples; none of our starburst galaxies had equivalent widths 
below log(£W(Ha)) < 1.51, while only two out of ten or- 
dinary galaxies had values greater than 1.55. 

This apparent disconnect between the starburst and 
WIM luminosities contrasts strongly with what is ob- 
served in ordinary star-forming galaxies. It is well-known 
that ordinarily, the WIM comprises about half of the 
total Ha luminosity in star-forming galaxies, regardless 
of galaxy lu minosity, morpho lo gy, or star- f ormat ion in- 
tensity (e.g. lWalterbosiri998( ). TOev et all (|2007[ ) stud- 
ied a much larger sample of star-forming galaxies than 
had previously been examined, confirmed this result, but 
with the notable exception of the starburst galaxies. The 



physical basis for the constant WIM fraction is not un- 
derstood, and the variation from this behavior by star- 
bursts may be an important clue. More work is needed 
to understand the radiation transfer of ionizing photons 
in both starbursts and ordinary star-forming galaxies. 

4. CONCLUSIONS 

Our multiwavelength comparison of starburst and non- 
starburst galaxies shows that on average, the SEDs of the 
starburst galaxies, defined by star-formation intensity, 
do not differ significantly from those of their ordinary 
star-forming counterparts. Within each band, individual 
starburst and non-starburst galaxies varied by up to 0.7 
dex in their flux ratios, and only in 24/zm and 160^m 
fluxes were the gaps between the mean flux ratios of the 
starburst and non-starburst galaxies larger than the sta- 
tistical uncertainty in each. As a result, our data do not 
support the hypotheses that ionizing radiation is either 
escaping or being absorbed by nearby dust in our sample 
of starbursts, in significantly larger fractions than would 
happen in ordinary star-forming galaxies. 

We do find evidence that dust in the starburst galaxies 
is somewhat hotter, as shown by the higher 24/70/im flux 
ratios and the extremely concentrated 24/im radial pro- 
file of starbursts. The coincidence of the Ha and 24^m 
radial profiles in starburst galaxies strongly suggests that 
the young stellar population traced by Ha flux is respon- 
sible for the increase in hot dust traced by 24/im in- 
frared flux. However, the overall distribution of the radi- 
ant energy emission in this Ha-selected sample shows the 
starbursts and non-starburst galaxies to be broadly simi- 
lar; while starburst galaxies appear to have substantially 
higher amounts of hot dust traced through 24/im flux, 
they have a correspondingly reduced amount of cold dust 
detected through 160/im flux, with total infrared fluxes 
comparable to their ordinary star-forming counterparts. 

We find no evidence that the low WIM fractions in 
starburst galaxies are caused by higher / esc . Further- 
more, there is no correlation of infrared excess with 
/wim j implying that dust absorption of the ionizing ra- 
diation also cannot explain this effect. Instead, an in- 
verse correlation of Ha equivalent width and /wim (Fig- 
ure [5]) strongly suggests that the WIM luminosity is inde- 
pendent of the starburst activity, and so the WIM con- 
tribution shrinks when the Ha luminosity is temporar- 
ily boosted by the starburst. This behavior contrasts 
strongly with the relationship between star formation 
and WIM luminosity in ordinary star-forming galaxies. 

In summary, this study is consistent with those that 
found low /esc from starbursts. Why these galaxies with 
extreme star formation have low / csc is not understood, 
nor do we know why the WIM luminosity correlates 
with SFR in ordinary galaxies, even though it is now 
apparent these do not correlate in starburst galaxies. 
These mysteries emphasize that our physical under- 
standing of radiation transfer for ionizing photons is 
extremely limited. We also caution that our sample is 
limited in parameter space, especially in terms of the 
magnitude of starburst activity, and additional studies, 
including those of more extreme systems, are needed to 
fully understand the circumstances that allow for high 
escape fractions of ionizing radiation. 
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